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ABSTRACT 


In a free electron laser (FEL) system, knowing the optical beam 
characteristics is of great importance. A beam may be comprised of higher-order 
modes due to the interaction with the electron beam, or from non-ideal 
operational conditions such as mirror distortions and misalignments, or from 
imperfect injection of the electron beam. 

In this thesis, the basic FEL theory is initially reviewed. The parabolic 
wave equation is then solved for the “fundamental” Gaussian mode and for 
higher-order modes. Working in rectangular coordinates, a complete and 
orthogonal set of solutions involving Hermite polynomials is found. When the 
wave equation is solved in cylindrical coordinates, we arrive at a set of solutions 
that contain Laguerre polynomials. The so-called Laguerre-Gaussian modes are 
analyzed. The evolution of these laser modes is also explored, yielding quite 
unexpected results due to their phase structure and orbital angular momentum of 
light. Lastly, we study a common case where higher-order optical modes appear, 
in order to quantify the tolerances of an FEL. 
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I. INTRODUCTION 


The Free Electron Laser (FEL) is a source of powerful, coherent radiation 
with unique attributes. Continuous tunability, high power, good efficiency and 
reliability make it suitable for a wide span of applications. The option of high 
power derives in part from the absence of a medium vulnerable to damage. 

Powerful coherent radiation was first produced when the microwave tube, 
which uses a beam of free non-relativistic electrons inside a closed microwave 
cavity, was invented in the 1930’s. Following this invention, the introduction of the 
open resonator enabled many steps towards the development of more powerful 
and efficient radiation sources. 

Based on its special attributes, the FEL is suitable for many diverse 
applications. It can be used in biology to study cell and molecular structure, and it 
can be used in the military, due to its high power output and its capability to 
operate in a broad range of wavelengths. A critical requirement for an FEL to 
have the desired results is good optical beam quality. We need to know the beam 
shape that can be produced under various operational conditions, so as to 
evaluate its performance. 

In this thesis, we analyze various high-order modes from an FEL, 
emphasizing cylindrical coordinates, which are motivated by the symmetry of the 
system. We present a brief review of a typical FEL system in the first chapter, as 
well as the relevant operating configurations. The next chapter introduces basic 
FEL theory, where we derive the pendulum equation, which describes the 
electron motion, and the parabolic wave equation, which describes the laser light 
evolution, from the classical electromagnetic theory. 

In the chapters that follow, we focus our attention on the optical beam. 
First, the basic optical theory is reviewed and the fundamental-mode solution to 
the wave equation is derived and analyzed. Although this solution is a good 
approximation for the output of many FELs, the actual operating parameters of a 
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system can generate laser light with higher-order modes. To analyze this effect, 
we initially review work with rectangular coordinates, and find the Hermite- 
Gaussian modes. Following this, new work is presented with cylindrical 
coordinates, which are more relevant to an FEL system, and we arrive at the 
Laguerre-Gaussian solutions. Analyzing these high-order modes, we find that the 
orbital angular momentum of light plays important role in the evolution of such a 
beam. 

In the last chapter, we present the results of a new set of simulations 
designed to investigate the tolerance of the Jefferson Lab FEL (Figure 1) to 
electron beam tilt. From these results, we see that higher-order optical modes do 
indeed appear as a consequence of the non-ideal operating conditions. 


Electron Gun High voltage 



Figure 1. Schematic of Jefferson Lab free electron laser After [1]. 
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II. FREE ELECTRON LASER SYSTEM 


In the Free Electron Laser, a beam of relativistic electrons, produced by 
an injector and accelerated by a linear electron accelerator, passes through a 
transverse, periodic magnetic field, called the undulator, oscillates and therefore 
radiates. Radiation can be captured in a cavity and used to induce new electrons 
to emit even more radiation. Alternatively, radiation from an external source can 
be amplified over a single pass, and even without an external source it is 
possible to produce radiation through the process of spontaneous emission. 
Therefore, it is informative to look at the possible basic setups, which can be 
used in order for a free electron laser system to function. 

A. FEL CONFIGURATIONS 

The configuration used in the first FEL by John Madey in 1977 is the 
oscillator, employing an optical cavity in a very similar way as in the conventional 
chemical or solid state lasers. Radiation is stored and amplified in the optical 
cavity, bracketed by one fully and one partially reflective mirror, the latter being 
used to outcouple the laser light. Coherence is established after many passes of 
the light. Figure 2 is a typical FEL oscillator, where the emphasis is placed on the 
resonator. 



electron source 
(gun) 


Figure 2. Oscillator Free Electron Laser From [2]. 
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In the amplifier configuration there is no optical resonator and it can be 
considered the simplest type of FEL. Light from a seed laser is amplified by the 
interaction with the electron beam, passing only once through the periodic 
magnetic field. Coherence in this case is set up by the seed laser. 

A special type of amplifier is the SASE FEL (Self-Amplified Spontaneous 
Emission) where there is no external source. Instead, the spontaneous radiation 
in the first part of the magnetic field is used as the seed laser. At present. It 
appears to be the most promising concept in the race towards achieving much 
shorter wavelengths (LIV,X-rays), since there are no seed lasers in this regime. 
Figure 3 shows the layout of a SASE FEL, where there are neither resonator 
mirrors nor seed laser, but where a significantly longer undulator is required. 



Figure 3. SASE Free Electron Laser From [2]. 

B. FEL COMPONENTS 

The major components of an FEL system are the electron beam source, 
the accelerator, the undulator, the resonator (for the oscillator) and the beam 
dump. There are many other elements that make an FEL system operational, like 
bending and control magnets to steer the electron beam in the desired path and 
several detecting devices to continuously monitor all the critical points, in order to 
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keep track of the operational parameters and to enhance the performance. Last 
but not least, an optical system must be utilized to transport the laser light in the 
desired area. 

1. Electron Beam Injector 

The electron beam injector is the original source of electrons in the FEL, 
and it plays a crucial roll in the performance of the system. It is usually a 
photoinjector, where a drive laser excites the electrons in a cathode, which is 
positioned inside a Radio Frequency (RF) cavity, and a high voltage is applied so 
as to accelerate the released electrons, allowing them to gain a modest amount 
of energy before entering the main accelerator. Instead of taking advantage of 
the photoelectric effect, much the same release of electrons can be obtained by 
heating the cathode. The electron beam that is produced from an injector, like the 
one shown in Figure 4, usually comes out in pulses at regular intervals. 



Figure 4. Electron Injector After [3]. 

2. Accelerator 

The electrons emerge from the injector with energy of about 5-10 MeV. In 

order to achieve infrared radiation from the FEL, we need a relativistic electron 

beam of about 100 MeV. To attain this energy level, the electron beam is passed 

through a Radio Frequency Linear Accelerator (RF LINAC). Superconducting 
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technology may be used as in the JLab FEL LINAC, shown in Figure 5. The 
accelerator uses an alternating electric field, which interacts with each electron 
bunch that passes through the RF cavity. For this reason, the electric field has to 
be fully synchronized with the electron bunches, so that it continuously 
accelerates the electrons. Also, the accelerator must have several stages of RF 
modules in order for the electrons to reach the required energy. 

The RF field inside the linear accelerator is generated by a klystron. To 
reduce the input power required, a special technique has been introduced. The 
electron beam that has traveled through the whole system can be driven back 
through the accelerator, so as to give most of its remaining energy, back to the 
RF field. To achieve this, we have to carefully adjust the timing of the beam, so 
as to be out of phase with the electric field of the RF cavities. This energy 
recovery procedure through the recirculation of the electron beam increases the 
overall efficiency of the system and additionally reduces the beam dump size. 



(c) 

Figure 5. JLab Superconducting LINAC (a), (b) A five-cell cavity, photo and 
sketch, (c) Part of a cryomodule which contains six five-cell cavities 
From [4]. 
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3. 


Undulator 


The undulator is the component of the FEL system that forces the electron 
beam to wiggle and radiate. It is composed of a group of magnets, as can be 
seen in Figures 2 and 3, arranged so that they form an alternating magnetic field 
along the axis of the undulator. The electrons leaving the accelerator with 
relativistic speeds, enter the periodic magnetic field of the undulator, where the 
Lorentz force drives them to oscillate in a transverse direction. The radiation that 
is emitted interacts with the light already present in the undulator, through the 
process of the stimulated emission. This interaction results in the amplification of 
the stored light. 

4. Resonator 

In the oscillator configuration, as the one shown in Figure 2, the FEL 
system makes use of an optical cavity with the undulator enclosed inside. This 
way, light is stored and amplified through many passes. Two mirrors bound the 
cavity: one of them totally reflects the light, while the other allows a percentage to 
come out. The light that emerges through this out-coupling mirror is the desired, 
usable laser light of the system. It is obvious that for maximum coupling 
efficiency between the electrons and the light within the resonator, the optical 
beam should be synchronized with the electron bunch. 

5. Beam Dump 

The beam dump is simply a notched piece of metal that absorbs the 
energy of the electrons at the end of their path, either just after the undulator or in 
the energy recovery configuration, after the linear accelerator. Even though 
electrons give energy to the optical field in order for the FEL to emit light, and 
probably back to the RF field, they still have energy to be removed. Therefore, 
the beam dump can heat up at megawatt rates and needs to be water-cooled. 
Furthermore, it is shielded so that no radiation can exit. A schematic diagram of a 
beam dump is shown in Figure 6. 
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Electron 


Figure 6. 
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Concrete 


Schematic diagram of the electron beam dump After [5]. 
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III. BASIC FREE ELECTRON LASER THEORY 


The basic idea behind the operation of a free electron laser is the 
extraction of energy from relativistic electrons, forcing them to radiate coherently. 
This energy conversion takes place inside the undulator, where the electron 
beam makes a wiggling, transverse motion, emitting spontaneous radiation in the 
forward direction. We can couple this motion to the electric field of the radiation 
to produce stimulated emission, thus amplifying the optical mode. We will now 
look at the mathematics that describe the attributes of a free electron laser and 
explain its operation. 

A. FEL PENDULUM EQUATION 

In order to understand the physics behind the electron motion and the 
subsequent radiation, we will first study their behavior inside the magnetic field of 
the undulator. For that reason, we will consider relativistic electrons entering a 
helical undulator field along the z-axis in the presence of optical field. As the 
electrons travel along the undulator, forces from both the static magnetic field of 
the undulator, and the oscillating electric and magnetic fields of the optical mode 
affect them [6]. 

The magnetic field of the undulator can be expressed in rectangular 
coordinates as: 

Bjj = B (cos (k^z), sin {k„z), o), (111-1) 

where B is the magnetic field strength, is the undulator wavenumber 

and is the undulator wavelength, being the distance of one complete cycle of 
the magnetic field in the undulator. 

The corresponding electric and magnetic field of the optical wave are 
given by: 
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E = E (cos y/, - sin y/, O), 

(III.2) 

5 = E (sin y/, cos y/,0P 

(III.3) 


with V^-kz + E the field amplitude of the optical wave in cgs units, 
k = (olc the optical wavenumber, co the frequency and (p the optical phase. 

Electrons entering these fields encounter a Lorentz force that, in the case 
of relativistic particles, is given by the equations 


d[yp) 

dt 


e 

me 


[e+Pxb), 


(111-4) 


^=-—P-E , (III.5) 

dt me 


where y= is the Lorentz factor, p= v/c is the dimensionless 

electron velocity, e = \e\ is the electron charge magnitude and m is the electron 

rest mass. Equation (III.5) is also known as the energy equation because it 
describes how electrons exchange energy with the optical field, and the electron 
energy is s=ymc^. 

We can now substitute the fields from (lll.1)-(lll.3) into (III.4) to obtain the 
spatial components of the force on the electron 


d[rh) 


dt 


= r£'(l-y^J(cos^,-sin^,0) + y5^5(-sin^^z,cos^^z,0)1, (III.6) 

me 


djrP,) 

dt 


e 

me 


E cos sin y/'yB sin k^z - cos ^„z)], 


(111-7) 
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where are the three components of the velocity and i+p^ j is 

the transverse velocity. Since we are dealing with relativistic electrons, we can 
take P^~\. Then, all terms of £'(l-/?^) are small compared to p^B and can be 

dropped. Integrating (III.6) and assuming perfect injection into helical orbits, so 
as to take the constants of integration to be zero, we obtain the velocity of the 
transverse motion 



Inymc^ 


(cos^^z,sin^^z,0). 


(III.8) 


We can write (III.8) in a more compact form if we introduce the dimensionless 
undulator parameter 


K= 


eB X 

rms o 


(III.9) 


with the rms value of the undulator magnetic field. Then 


=-—(cos^^z,sin^^z,0). (111.10) 

r 

In (III.10) note that for the helical undulator, = B and for a linearly polarized 
undulator B^^^ =b/>/2 . 

At this point, we are able to look at the basic transverse motion of the 
electrons, by making the assumption that they travel roughly with the speed of 
light along the undulator. Then , and taking the electron velocity to be 

constant we can write p^ct. Using a>^ = p^ k^c as the electron 

oscillation frequency, (111.10) becomes 

Pj^~ -y (cos( 0 /),sin(ft>„r),O), (111.11) 
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and integrating we get 



K /I 

— T^(-sin(®„r),cos(®„r),0) 

Y V / / 


( 111 . 12 ) 


with x_^[t) = x[t)i+y[t)j. This equation approximates the transverse motion of 

the electrons inside the undulator. From this expression and for typical systems, 
we find that electrons move in the transverse directions in a distance of about 
25 jum, the oscillation frequency is in the microwave regime, but due to the 

relativistic Doppler effect, the radiation frequency is multiplied by a factor of . 
resulting in lasing at much higher infrared frequencies. 

Next, we want to study more carefully the microscopic motion of the 
electrons, without making the approximations of the previous paragraph. For this 
reason, we will use the transverse motion we found in equation (III.10), and plug 
it into the energy equation (III.5) to get 

Y = ^=-—e( 5^cost//-^ sini//\ (111.13) 

Introducing the electron phase ^={k + k^^z- cot and recalling that 
W = kz-(ot-t-(p^ (111.13) is rewritten as 

/= ^^^cos(i^+^). (111.14) 

yme 

The last equation describes energy exchange, because it tells us whether 
electrons lose or gain energy, based on the sign of the cos{^+(p) term. For 

instance, when cos(i^+^)<0 then y<0, thus electron energy 

decreases, which is the desired condition for the FEL to operate, so that 
electrons give energy to the optical field. 
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Electron phase ^ is a way to relate electron position to the co-propagating 
optical field; it corresponds to the position of the electron in a section of the 
electron beam one optical wavelength long. Therefore, it is reasonable to 
introduce the electron phase velocity, defined in the following way so as to be 
dimensionless 




(111.15) 


where L = NX^ is the length of the undulator, N is the number of undulator 

periods and the symbol (..) indicates the derivative with respect to the 
dimensionless time t , defined by z = ctlL. This normalizes time to the length of 
the undulator, so that r = 0 ^ 1 from the beginning to the end of the undulator. 

The dimensionless electron phase velocity v is also known as the 
“resonance parameter”. For an electron that satisfies the resonance condition 
v=0, the electrons and optical field have optimum energy exchange. In this 
case, c/?^ = kc/(k + k^y Assuming again relativistic electrons, we have «1 and 
k:s>k^. Therefore, X= X^ and laser frequency where a>^ ~ Inc/X^ is the 

electron oscillation frequency. Then, the laser wavelength is 

A 


From (111.10), we find that 0[ = and so 





(III.17) 


For relativistic electrons with / »!, equation (111.17) can be approximated by 
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(111.18) 


A«1- 


\ + K^ 


Combining equations (III.16) and (111.18) we find the following expression of the 
resonance condition for the laser wavelength 


A = 



(111.19) 


This resonance relationship makes obvious that an FEL can be continuously 
tuned by changing physical parameters, like the electron energy £ = /mc^ or the 
design of the undulator. Since the undulator parameter is typically «1, it is 
clear that optical wavelength is much smaller than that of the undulator by a 
factor of approximately . 

We are able now to continue our analysis, so as to arrive to a compact 
expression for the electron microscopic motion, taking first the time derivatives of 
the electron phase velocity v, defined in (111.15) and that of the longitudinal 
velocity from (111.18) 

v^L(k + K)i},, (111.20) 
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Then, we combine these results to find v = A7rNYlY and finally substituting f 
from (111.14) we get 


V 



ATiNeKLE 

2 2 

Y me 


cos 




( 111 . 21 ) 


It is now reasonable to collect terms and define \a\ = 4;rNeKLEjY^mc^ as the 

dimensionless laser field amplitude, so that (111.21) becomes 
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( 111 . 22 ) 


0 00 

v=<^= |a|cos 


{C+(p)^ 


which is the FEL “pendulum equation” that describes the microscopic motion of 
each electron in the reference frame of one optical wavelength. In this reference 
frame, the motion of the electrons can be taken to be similar to a classical 
pendulum. The dimensionless laser field amplitude |a| is a measure of the 
strength of the optical field and expresses the rate of the electron bunching. 
When I a| » the laser has strong fields and it is near saturation, while for |a| < tt 
the optical field is weak and there is no overbunching of electrons. 

Equation (111.22) implies that electrons inside the magnetic field of the 
undulator move in an equivalent way to the classical pendulum. Phase space 
plots are commonly utilized to visualize the motion of a pendulum, thus they can 
also be used to describe the motion of the electrons within an optical wavelength, 
showing the evolution of electron phase and phase velocity (analogous to 
position and angular velocity for the classical pendulum). Since phase velocity is 
related to the energy, phase space plots give an informative picture of the energy 
exchange between electrons and the optical field. 

The following figures show the phase space paths for 20 sample 
electrons, injected with the same initial phase velocity, and having uniformly 
distributed initial phase. The FEL phase space paths can be found from the 
expression [6] 

v" = vt + 2|a|[sin(4'+ ^)- sin(4; + (p)\ , (111.23) 

while the separatrix is given by 

= 2|a|[l + sin(4'+^(9)] i.e., v= ±-^2|a| [l + sin(4'+^)]^, (111.24) 
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which defines a closed contour inside of which are the closed orbits, and outside 
of which are the open orbits. It passes through the unstable fixed points 
v)= (-;t/2,0) and {^,v)={37rl2,Q), and its peak to peak height is easily seen 

to be 4|a|'^^. An electron moving in an open orbit corresponds to a pendulum that 

swings over the top, whereas an electron that follows a closed orbit corresponds 
to a pendulum that swings back and forth. 

In Figure 7 it is plotted the evolution of 20 sample electrons starting at 
resonance ( =0), during the time period from r = 0 (appearing as yellow dots) 
to r = l (appearing as red dots). All paths are inside the separatrix, thus all 
electrons follow closed paths. The symmetry that appears, indicates that half of 
the electrons gain energy from the optical field and their paths shift up, while the 
other half lose energy. This is represented in the same figure from the graph of 
the net gain G(t), which is almost zero. Gain is defined as 

where p{t) is the optical power at time r and is the initial power. The optical 
phase evolution (pif) is plotted in the same figure. 
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Figure 7. Phase Space evolution of Electrons at Resonance (v^=0), 
corresponding optical field Gain and Phase evolution. 

Contrary to the nearly zero net energy transfer when electrons are injected 
at resonance, we can arrange to have electrons give net energy to the optical 
field. Injecting electrons slightly off resonance, as shown in Figure 8, we notice 
that while some of them perform open orbits, at the end of the undulator when 
r = 1, most electrons have moved to a lower value of phase velocity, resulting in 
a net loss of energy. In this case, a very important effect takes place near (^= tt , 
namely “bunching” of the electrons and a final gain of more than 10%. 
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Figure 8. Phase Space evolution of Electrons injected off resonance 
(= 2.6), corresponding optical field Gain and Phase evolution. 
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B. DIFFRACTION OF LASER BEAMS 


We will now study the propagation of laser light through a sourceless 
medium, in order to be able to understand and describe how light travels and 
moreover, how is affected by diffraction. We start with Maxwell’s equations, 
assuming for our case an isotropic and homogeneous medium. When there are 
no sources, that is to say the charge density p = Q and the current density 7 = 0 , 
we write the four Maxwell’s equations in cgs units as 


V£ = 0, 

(111.26) 

¥•5 = 0, 

(111.27) 

Vx5 =-, 

c dt 

(111.28) 

Vx5 =-, 

c dt 

(111.29) 


where c is the speed of light. Taking the curl of both sides of equation (111.28) we 
have 


Vx 



r 


= Vx 


V 


1 dB^ 

c dt ^ 



Applying a vector identity on the left-hand side and substituting VxB from 
(111.29), we get 


v(v-£)-V"£ = --— 

^ ^ c dt 


^\dE^ 

yC dt j 


(111.30) 


Then, since V£’ = 0 from (111.26), we get 
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f 


V" 


c" de ^ 


E = 0, 


(111.31) 


which is the wave equation that describes light propagation. Electric field E = Ee, 

where s is the polarization vector, is a function of position and time, thus 
equation (III.31) implies four second-order derivatives in (x,r). Due to the 

coherence of the laser light, electric field E is slowly varying in the direction of 
propagation over an optical wavelength. Likewise, E is slowly varying in time 
over the optical period. Then, we can write the magnitude of the complex electric 
field in the form 

E (x, r) = cr (x, r) 

= (111.32) 


where kz-cot is the carrier wave, a = \a\e‘‘^ is the complex optical field and 
x={x,y,z). Next, we want to substitute the electric field (111.32) into the wave 
equation (111.31). After taking the required derivatives, we find 


V" 


oc +1 


( 


dz^ 


a ^da 
+ 2ik- 


^ 1 


dz 


k a \ - 




f d^a 


de 


2ico 


da 

dt 


co^a 


i{kz-cot) 


= 0 , 


1.33) 


with = d^/dx^ +d^ldy^ the transverse Laplacian. We can simplify this equation 
by dropping out the second derivatives, because laser field amplitude |a(x,r)| 
and phase (p{x,t) are slowly varying. This means that 

,1 I II ^ 1 /111 o/i\ 

—^ = k\a\ , — — co\a\ , —= cocp and —= kcp. (111.34) 

dz dt dt dz 
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From d\a\ldz = k\a\ we find that d^\a\ldz = kd\a\ldz and similarly from 

d\a\ldt= co\a\ we find that d^\a\ldt^= cod\a\ldt, SO that we are justified in 

dropping the second derivatives in favor of the first derivatives in (111.33). In the 
same equation, the last term in each parenthesis cancels out, since 

k = CO ! c^k^ = 0 / / k^a = co^a/c^ . Then, after multiplying by 
equation (111.33) becomes 



^ d 15^ 

V" +2ik 

- 1 - 

± 1 

\dz c dt J_ 


cr(x,r) = 0. 


(111.35) 


In order to simplify the last expression, we introduce a new variable 
u = z-ct and again use the dimensionless time T = ctlL (where L is the 

propagation range). We can now rewrite the operator d/dz + d/dt that appears 
in (111.35) taking the necessary partial derivatives 

d du d dr d 

dz dz du dz dr > (111.36) 

\ d \ du d ^ \ dr d 
c dt c dt du c dt dr 


and since 5a/5z = l, 5r/5z = 0, duldt = -c and 5r/5r = c/L, the operator 
dldz + c~^ d/dt is just U^d/dz and consequently the wave equation (111.35) is 
now written 




V" + 


7.ik d ^ 
L dzj 


a(x,r) = 0, 


(111.37) 


which is called the “parabolic” or “paraxial” wave equation, because we have 
made use of the paraxial approximations in (111.34). 
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At this point, it is convenient to use dimensionless transverse 
coordinates (x,y) along with the normalized time r. If we rearrange the terms in 
equation (111.37), we can rewrite it as 


l^J 


V" + 


d 

dr 


ct'(x,r) = 0. 


(111.38) 


where we notice from the first term, which is the diffraction term of the wave 
equation, that the transverse extent of the laser beam, together with the 
wavenumber k = 2nlX and the propagation range L, affect the diffraction of the 
beam. The equation takes a simpler form if we normalize the transverse 
coordinates using the relationships 


X = 


f 1. \ 


v2Ly 



y = y 


yiL) 


(111.39) 


In (111.39), the normalization factor is related to the characteristic mode 

waste radius by 


= 




(111.40) 


with Z^ = L in our case, where is the Rayleigh length, defined to be the 
propagation distance in which the beam doubles its area and it is determined by 
the curvature of the mirrors. Finally, using only dimensionless coordinates and 
utilizing the dimensionless transverse Laplacian, we arrive at the dimensionless 
parabolic wave equation, written without the tildes for convenience 




4 dr j 


a(x,r) = 0. 


(111.41) 
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From this last wave equation, we notice that for a beam with transverse 
dimensions of the order of unity or less, we have to take into account the 
diffraction, since it affects the beam a great deal, over the range r = 0 ^ 1 , while 
for a much wider beam we can ignore it. Diffraction of the wavefront causes both 
the amplitude and the phase of the optical field to evolve. 


C. THE FEL WAVE EQUATION 


We have found, so far, a mathematical description for the laser light that 
propagates in absence of sources and also for the electron motion in the 
presence of light. What is left, in order to generalize, is to describe the optical 
wave evolution, taking into account the interaction with the current source due to 
the electron beam. We start with the full wave equation [6] 


f 


1 a 


2 A 


- 


A{x,t) = -^J^{x,t) 


( 111 . 42 ) 


where A[x,t) is the optical vector potential and J±{x,t) is the transverse current 

density due to oscillations of the electrons passing through the undulator. Then 
we are able to find the optical electric and magnetic fields from the vector 
potential using the relationships 

— 1a4 — — 

E = -, 5 = VxA. ( 111 . 43 ) 

c dt 

Once again, we will make use of the fact that laser light is coherent. Thus, we 
take A{x,t) to have a slowly varying envelope in the direction of propagation. In 

addition, we take the electric field to vary slowly over a laser wavelength and to 
vary slowly in time, relative to the optical frequency. Assuming a complex laser 
electric field E = \E\e'‘^ , the vector potential can be written [6] 
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(111.44) 




where kz - cot is the carrier wave and s is the laser polarization vector, which is 
€ ={-i,\,0) for the case of a helical undulator producing circular polarized laser 

light. If we assume again slowly varying wave amplitude and phase and use 

(111.44), the wave equation (111.42) becomes 


s e 


[kz-cot) 


V\+2ik 


(d 15 ^ 

- 1 - 

V c dt j 


Atc — 
E = — -J± 
c 



r_ . . 

( d 

1 


I 

< 

+ 

— + — 


ydz 

c dt J 




(111.45) 


We will now use the same technique as we used in (111.36), introducing the new 
variable u = z-ct and r = cr/L, where L is the length of the undulator, so as to 

rewrite the operator d/dz + c”' d/dt. Then (111.45) becomes 




vi + 


2ik 5 ^ „ Aztk — 

- E =- J L -s 

L dr) c 


e 


-i{kz-cot) 


(111.46) 


This last equation is the parabolic wave equation with a source current 
7±, due to the bunched electrons. The source current, which is due to all of the 
electrons, can be expressed as 


where is the three-dimensional Dirac delta function, defined as 

^(^)(r) = ^(x)^(y)<J(z) and r.{t) is the position of the ith electron at time t. The 

transverse velocity found in (111.10) can be expressed as 
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(111.48) 



(cos k^z, sin k^z, O) = iRe 



Plugging (111.48) into the source current relationship in (111.47) and then the 
resulting source current into the parabolic wave equation (111.46), we have 


L dv 


E = -AmeKk p{x,t) {e . 


(111.49) 


In this equation, p{x,t) is the electron beam particle density in a small volume 
element dV 


p{x,t) = ^(0) ’ 


(111.50) 


and with the symbol (.) we denote the average over sample electrons. 

Employing again the dimensionless time r = cr/L, where L is the length of the 

undulator, and defining the dimensionless complex laser field a = |a|e'^, with the 
laser field amplitude |a| = ATzNeKLEjylmc^ , as used in (111.22), the wave equation 
becomes 




-itvi+A'' 

2k dv j 


{x,t) = - 



(111.51) 


where we have introduced another important PEL parameter, the dimensionless 
current density j = S7r^Ne^K^L^p/ylmc^ , and we have made the reasonable 
assumption that y ~y„ during the whole interaction. We want now to make the 
wave equation completely dimensionless, thus we bring in again the 
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dimensionless transverse coordinates x = x(^/2L)^, y = y[kf2LY^ and the 

dimensionless Laplacian operator. Dropping the tildes, the FEL wave equation is 
finally written 




--V"+—^ 
4 dr 


1.52) 


{x,t) = - 



In the case when the diffraction can be neglected, the last equation is 
simplified to 




(111.53) 


From this expression, we realize that the electron beam should be bunched, 
since the factor is a measure of the electron bunching. For instance, if the 

electrons are bunched at around the laser field amplitude goes up, 

whereas if the bunching occurs at 1 ^= 0, the field amplitude goes down. Equation 
(1111.53) also suggests that the laser field is related to the dimensionless current 
density j , so that when j <7r \Ne are in the low gain regime, while for j the 
coupling between the laser light and the electron beam is large and we have high 
gain. 


D. GAIN 

In order to better understand and describe the energy exchange between 
the electron and the optical beam, it is handy to make use of the gain we defined 
in (111.25), which is the fractional change of the optical field power during a single 
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pass. In terms of the optical field amplitude, gain can be written as 
G(r)= ' ^ ^ , ° =-^—1 or G(r)«2| ^^-11 , (111.54) 

\ J 

for low gain and small change in |a(r)| from the initial field amplitude a^. 

It is particularly useful to study the gain in the weak optical field state, 
because it gives us information about the performance of the FEL. For that 
reason, we will derive an expression for the gain, assuming weak optical fields, 
TT and low gain, so as to take |a|»a^and for the whole interaction. 

Applying the perturbation theory, we start by expanding the electron phase 
and the electron phase velocity v in powers of the initial small field amplitude 

+ + + v= +++ (111.55) 

where is the zeroth order in a^, is the first order, and so on. We can now 
plug (111.55) into the pendulum equation (111.22) to get 

H—j. (111.56) 

00 o 

From the zeroth order in a^ we have which gives the solutions 

= (^^ + v^T = 0 and . From the first order in a^, using the zeroth order 

solutions and integrating over time t , we find the solutions 

= _-^[cos(C + y,r)- cos(c)+ v„rsin(c)], 

J == ^[sin(4; + v^t)- sm(4;)]. (111.57) 
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We can now make use of these results in the similarly expanded wave equation 
(111.53) 

a = - j(e-‘^)= - . (111.58) 

Taking into account only the zeroth and the first order, we can approximate the 
last equation 

a « - , (111.59) 

but since = J = 0, (111.59) is written 

a « . (111.60) 


We can plug in from (111.57) and integrate, recalling that the complex optical 
field is a = lale'^, to find 


|a(r)| 


( ^2-2cos(v„r)-V rsin(v„r)')^ 

- ^2..- jj""' 


(P(.T^) = j 


2sm{yj)- vj(\ + cos(v^r))^ 
2v' 


+ • 


(111.61) 


Now, using equation (III.54), we are able to arrive at the low gain equation along 
the undulator [6] 


( 2-2cos(v r)-v„rsin(v r)^ 
G(r) = ;|- ^”^3 


1 I 

J- 


(111.62) 
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This equation indicates that the gain, in the weak field approximation, is 
proportional to the dimensionless current j and the initial electron phase velocity 
v^. Plotting the final gain spectrum at r = l in a single pass through the 
undulator, as shown in Figure 9, we can verify the results we found in the phase 
space plots of Figures 7 and 8. Particularly, when the electrons are ejected 
exactly at resonance v„ = 0, then there is zero net gain, meaning that the 
electrons lose and gain the same amount of energy during their interaction with 
the optical field inside the undulator. We also note, that the gain curve is anti¬ 
symmetric in . The peak gain of about 13% is obtained when the FEL operates 
off resonance at «2.6, while the peak absorption occurs at «-2.6. These 
results are consistent with the phase space analysis, suggesting that we must 
operate the FEL slightly off resonance in order to have the best exchange of 
energy from the electron beam to the optical mode and hence the maximum 
gain. 



Figure 9. Single pass gain versus initial electron phase velocity at r = 1. 
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IV. OPTICAL THEORY 


Up to this point we have examined the interaction of the electron beam and 
the optical field inside the undulator, arriving at some important expressions. We 
have also introduced the gain G(r) , a useful tool to measure the performance of 
this interaction. In this section, we will focus our study only on the optical field, 
propagating in free space. Equation (111.52) describes the optical field evolution, 
in the more general case of the presence of source. However, when light is 
traveling in a sourceless medium, specifically inside the undulator while there is 
no electron beam present and outside the undulator, this equation reduces to 
equation 





d ^ 
dr ^ 


a(x,r) = 0. 


(IV.1) 


In the following sections, we will find solutions to this wave equation, starting with 
the one that describes typical optical waves produced by a Free Electron Laser. 
We call this solution the “fundamental”, and it is the lowest-order mode. Higher- 
order modes, which we are going to analyze later, may also appear. 


A. FUNDAMENTAL SOLUTION TO THE WAVE EQUATION 

An ideal, realizable, laser beam cross sectional shape is Gaussian. We 
therefore seek a Gaussian solution to the wave equation. A trial solution of this 
form can be written 

a(r,r) = a^p(r)e(IV.2) 

where a(r,r) is the dimensionless optical field, and is the optical field 
amplitude at the point (x,y,r)= (0,0,0), +/ . The functions p{t) and q{T), 
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as yet unknown, will be determined by the wave equation (IV.1) and the initial 
conditions imposed by the desired Gaussian shape of the beam, 


a{r,0) = a^e , 


(IV.3) 


where = is the dimensionless Rayleigh length, and L is the propagation 
distance. These initial conditions give 

p(0) = l and qiO) = z,. (IV.4) 


In order to find the expression we are looking for, from the trial solution (IV.2), we 
have to find the functions p(t) and qir). Before putting (IV.2) into the wave 
equation, we switch to cylindrical coordinates, due to the symmetry of our 
problem. Then, the wave equation (IV.1) becomes 


i 1 d 

^ d^ 

4 r dr 

K drj 



a(r,r) = 0. 


(IV.5) 


Now we insert the trial solution (IV.2) into (IV.5). The first derivative of a(r,r) 
with respect to r is 


dr q(T) 


(IV.6) 


from which 





f 2 

r 


. dr) 

q{T) 

\qij) 


Finally, we calculate the time derivative of a(r,r) 


(IV.7) 
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f 


(IV.8) 


5a 

5r 


= a„p(r)e 


^(r) 


1 

'v p{t) dr 


+ 


2 7 A 

r dq 
q^ir) dr) 


We can now go back to the wave equation (IV.5) and employ these derivatives 


i 1 d 
4 r dr 


5a^ 5a ^ 
r— + — = 0 => 
drJ dr 


-la 


° q(T)^ Kq(T) 


q(T) 


- 11 + a^p(T)e 


1 dp 


+ • 


dq 


\p{T)dT q(T)dTj 


= 0 


(IV.9) 


Dividing through by a^p{T)e we get 


r 1 1 dp r dq 

-i - +1 -+- - + ^-^ = 0 

q{T) q{T) p{t) dr q (r) dr 


(IV.10) 


Multiplying by -iq^ir) and rearranging, we can write the last expression as 


dr 


1 i dp 


Kqir) p(T)dT) 


= 0 . 


(IV.11) 


In order for this equation to hold for all r at any time r , the quantities in the two 
parentheses must be identically zero, and therefore 

-l-/^ = 0, and 
dr 

_= 0 

q{T) p(t) dr 

From equation (IV.12), after integrating, we find q(T) = iT + q^ , and from the initial 
conditions (IV.4) we see that ^(0) = , so that 


(IV.12) 

(IV.13) 
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q{T) = iT + z„ . 


(IV.14) 


Using the last result in equation (IV. 13) we see that 

dp ^ ip(T) ^ p(t) ^ 

dr iT + z„ + 

J -}—dp = f ^ 

^ P{t) ^-t'+iz, 

ln/7(r) = -ln(-r'+«zJlo=-ln(-r + «zJ + ln(/zJ . (IV.15) 

and finally 


p{T) = 



(IV.16) 


We can now substitute p{t) and q{T) into (IV.2) and write down the solution to 
the parabolic wave equation as 


a(r,r) = a„ 


1 + ^V 


(IV.17) 


This solution can be written in a more informative way, expressed in terms of the 
dimensionless beam width w{t) , which is defined as 


w{t) = 




or 



(IV.18) 


We note that at r = 0, the beam width or waist in this case (the minimum value of 
the beam width), is just w(0) = . Of course, the beam waist need not be 
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at the beginning of propagation distance, thus a more general expression for the 
beam width at any point is 


w{t) = 





(IV.19) 


where is the dimensionless position of the beam waist (normalized, as usual, 
to the propagation distance L). A cross-section of a propagating Gaussian beam 
is shown in Figure 10. Physically, the beam width is the distance where the 
optical field drops to Ije of its maximum, as shown in Figure 11. 



Figure 10. Cross-section of a propagating Gaussian Beam. 

We are now ready to go back and modify the solution (IV. 17). To make the 
calculations easier we let u = \ + iT/z „. Then, equation (IV.17) becomes 


a(r,r) 



(IV.20) 
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Figure 11. Gaussian beam profile. 


and since the complex conjugate of u is w* = 1 - , we can write (IV.20) as 


i(r,r) = a„-f 


u 




hr^lzp 


L(r,r) = a„-e . 


(IV.21) 


The complex number u can be written in polar form as 


[ f /arctanj^^l 


(IV.22) 


and since uu* = 1 + {jlZof, equation (IV.21) becomes 
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a(r,r) = a„ 




^oTK^Aofl 

L J ^ L J 


Vi + (^A.T 


/ \ ^ 


(IV.23) 


We now insert the beam width w(r), defined in (IV. 18), into the last equation, 
which becomes 


/ \ W 2.x 2,2 -/arctan V 

a(r,r) = a„— 2 _ g g=:» g v/. 
^ ^ >v(r) 


(IV.24) 


where we also used . Finally, we can write the solution to the parabolic 

wave equation (IV.1) more compactly, realizing that the last two exponents in 
(IV.24) are pure imaginary. So, we let 


= ^ ~ arctan 


2 , 2 
Z„ +T 




v-^o y 


(IV.25) 


and write the fundamental solution as 


w(t) 




(IV.26) 


B. PROPERTIES AND PROPAGATION OF GAUSSIAN BEAMS 

The Gaussian beam is the lowest-order mode. A beam of this profile can 
be fully described by the beam waist and the position of the beam waist . 
We will study the behavior of the optical beams described by the solution (IV.26), 
in order to better understand the evolution of their properties as they propagate. 
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First, we note that, at the beginning of the propagation, when r = 0, the 
beam has the Gaussian shape that we chose with the initial conditions, namely 

a(r,0) = ay^». (IV.27) 

As time increases, the beam width defined in (IV. 18) gets bigger and so 
the transverse area of the beam expands. As a consequence, for a fixed distance 
from the center (fixed r), the field strength goes down as the optical field 
propagates. Looking at the center of the optical mode, where r = 0, the field 
decreases continuously, as it is expected, since a(0,z-) goes like ~ l/r. 

Gaussian modes have the very important feature that they do not change 
their profile as they propagate, except that the beam radius varies. In other 
words, a Gaussian beam stays Gaussian. This is the reason why we where 
looking for a solution to the wave equation of this form. 

Figure 12 represents a computer simulation of the propagation of a laser 
light in the fundamental mode. In this figure is plotted a cross-section of the 
Gaussian mode as it propagates from r = 0 to t = \, together with the transverse 
cross-section at the beginning and at the end of the propagation distance. 



Figure 12. Cross-section of a propagating Gaussian Beam with r^ = 0 and 

z„=0.3. 

Additionally, in Figure 13 is shown the surface plots of a Gaussian mode, at 
the beginning of the propagation when r = 0, at r = 0.5 and at the end of the 
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propagation distance when r=l (note the difference in the vertical scale 
between the three graphs). From all these plots, we are able to visualize what we 
analyzed previously, specifically the fact that the light beam expands as it 
propagates and also that it suffers attenuation. 




(a) 


(b) 



(c) 

Figure 13. Surface plots of a Gaussian Beam evolution at times (a) r = 0, 
(b) r = 0.5 and (c) r = 1. Scales in z-axis are different. 
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V. HIGHER ORDER MODES 


In the previous chapter, we considered only the fundamental Gaussian 
solution to the wave equation, which has the simplest intensity profile and is a 
good approximation to describe a laser beam. However, if we want to generalize 
for non-ideal beams and include the possible variations in the transverse beam 
profile, we have to find other solutions that correspond to higher order modes. 
These solutions to the parabolic wave equation have similar properties to the 
Gaussian. Furthermore, they form a complete and orthogonal set of functions, 
called the modes of propagation, so that every arbitrary optical beam can be 
expanded in terms of them [7]. We will find that the fundamental Gaussian beam 
is just the lowest-order mode of the generalized solution. 

As one can guess, the need for studying the optical beam more 
extensively arises from both experimental results and simulations. Generally 
speaking, whenever an FEL operates under non-ideal conditions, the laser beam 
may contain more than one higher-order modes. Such conditions could be 
vibrations of the mirrors in a portable system, electron beam shifts and tilts, 
mirror imperfections or deformations, and so on. For instance, vibrations can 
cause the electron beam to be injected slightly off-axis and to deform the mirrors, 
while heat can change the radius of curvature and shape of the mirrors. 


A. HERMITE-GAUSSIAN MODES 

The Hermite-Gaussian solutions are derived from the parabolic wave 
equation in rectangular coordinates, starting with the more general trial solution 
[ 8 ] 


a(r,r) = a„ g 


( 


X 

w { t ) 




w { t ) 


P{t)+ 


Zo q(z) 


(V.1) 
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where the functions p{t), q{T), ^(x/w)and /?(y/w)are found following a 

procedure similar to the one used for the fundamental mode. The same initial 
conditions as before require that 

/>(0) = 0, ^(0) = 1, (V.2) 

s{xlw„)=\, h{ylw^)=\. 

Using (V.2) and substituting the trial solution into the wave equation, we arrive at 
equations that have the form of the Hermite equation and can be solved in terms 
of the Hermite polynomials. After finding all the four unknown functions p, q, g 
and h, we write the m,n mode of the Hermite-Gaussian beam as 



H 

rviy] 

w 

\ J 

n 

w 

y J 


a (F,r) = a„^ H„ 

m,n \ ’ / o n 

w 

(V.3) 

where the phase is 

= - (OT + n + l)arctan 






2 2 
w 


K^oJ 


(V.4) 


and m,n are positive integers identifying the mode. In (V.3), H^and H„are the 
Hermite polynomials with m, n the non-negative integer indices. These indices 
determine the shape of the beam profile in the x and y direction respectively. The 
first ten of the Hermite polynomials are listed in Table 1. It is obvious that when 
m = 0 and n = 0, the Hermite-Gaussian beam is just the Gaussian beam, which 
is why it is called the fundamental mode. A few examples of higher order 
Hermite-Gaussian modes are shown in Figure 14, where we note that m nodes 
appear running vertically, while n nodes run horizontally. 
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//j(m)= 2u 

H^{u)=4\x"-2 

H^(u^= 8u^ - 12a 

//4(a)=16u"-48a'+12 

//5(a)=32u^-160a^ + 120a 

H^{u)= 64u® - 480a^ + 720a" -120 

Hj (a)= 128u" - 1344a^ + 3360a" - 1680a 

H, (a)= 256u' - 3584a® + 13440a" - 13440a" + 1680 

H, (a)= 512u‘' - 9216a" + 48384a" - 80640a" + 30240a 

Table 1. Hermite polynomials 



m = l, a = 0 m = 0, n = 3 m = l, n = l 



Figure 14. Hermite-Gaussian modes in the transverse plane x, y. Red 
indicates highest intensity, while blue indicates least intensity. 
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B. LAGUERRE-GAUSSIAN MODES 


The optical cavity of an FEL, utilizing circular mirrors and a round, aligned 
electron beam, has cylindrical symmetry. This fact suggests that we should use 
cylindrical coordinates to solve the dimensionless parabolic wave equation (IV.1), 
which is now written 


d^a 1 5a ^ 1 d^a ^ • 5a 
5r^ r dr d&^ dr 


(V.5) 


since the dimensionless transverse Laplacian in cylindrical coordinates is 


5 ^ 1_5 1 5 " 

5^2 ^ ^2 


(V.6) 


Motivated from the success in finding the Hermite-Gaussian modes trying the 
solution (V.1), we will now put a similar exponential term in the trial solution for 

9 

the wave equation, adding e to match the single valued periodic boundary 
conditions, where I is an integer. We also guess a power series in r , and since 
we anticipate diffraction that suggests the use of r/w dependence, we try the 
solution 


a(r,6»,r) = a^ 


V 2 r 
w(r) 




^ 2r^ ^ 


w (r) 


Zo ‘3(z) 


(V.7) 


where , P(r) and ^(r) are the unknown functions that we will try to 

find next. Although we use cylindrical coordinates, we still require that at r = 0 


a(r,i9, 0 ) — 


(V.8) 


42 



which gives i[i/lw") = 1, P(0) = 0 and ^(0) = 1. In order to put the trial solution 
(V.7) into the wave equation (V.5), we first need to calculate the derivatives of the 
function liir/w") using the chain rule 


diilrlw) di{2r/w-) 4 ^ 

ar a(2rVw^) ar 

where we used the prime to represent the first derivative of i with respect to its 
argument . Similarly, using the double prime for the second derivative, 

we find 


dHdrIw) 

ar' 


a 

f4r ^ 

_ 4 

(m 5<”'i 

i -\-r - 

dr 

1 w J 


1 ar J 


4/^ . 16r^ 


The time derivative of £[irl\v) is 


(V.10) 


d£(2r"jw") _ di[l/lw") a(2rVw^) _ p,/ 

dr dilr^jw^^ dr ^ ^ dr 


(V.11) 


with (..) indicating the derivative with respect to the dimensionless time t . Now 
we are ready to calculate the derivatives of the trial solution (V.7) 
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(V.12) 


(V.13) 
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(V.14) 
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(V.15) 


(V.16) 


Next, we want to substitute the trial solution (V.7) into the wave equation (V.5). 

/ /-/ x/ - P+—-(W 

Bringing all the derivatives together and dividing through by a(^2/w)e^ 
we have 


5^a 15a 1 5^a , 5a 

- 7 +-+ -r- 7 + 4«- 

5r^ r dr 63^ dr 
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r ir q 
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Zo q Z„q 


only dependence 


^-^-/P = 0. (V.17) 

Z,q w 

^ -V-^ 

no r dependence 


We want this equation to have solutions for each transverse plane, or in other 
words, for each fixed time t . Thus, we need to set the r and £ dependent terms 

equal to a constant, while the dependent terms should be equal to zero 
because equation (V.17) must be valid for any transverse position. Eventually, 
the terms with no r dependence must be equal and opposite to the constant 
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value of the terms dependent on r and i, in order for the equality to be true. 
Hence, we have to solve the following three equations simultaneously 


4r^ e 2 e Aiwr^ l' 
w'^ I I w^Zg q £ £ 


2 • 2 

r ir q 

2 2 +-^ 

Zo q Z„q 


= 0 , 


1 + 1 
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-iP = ^, 


(V.18) 

(V.19) 

(V.20) 


where we used - as the constant for the r and £ dependent terms. We will 
start from (V.19), since it is straightforward 


2 • 2 

r i r q 

+7 
Zo q Z„q 


= 0 


( 1 

Zoq^^z, 


= 0 


1 ° ^ 
- \-i q = 0, 


(V.21) 


where for the last step we use the same argument as we did for the derivation of 
the fundamental mode, namely that the equation must hold for all r. After 
integrating, and since ^(0) = 1 from the initial conditions, we find 


q(T) = — + l. (V.22) 

From the definition of the beam width w(t) = in (IV.18), we have 

0 

w = r/wz^. Rearranging terms in equation (V.18) we write 
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(V.23) 
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where the terms in the parenthesis, after substituting q{T) from (V.22) and using 
= w/, become 


-?-^- 7 "^-~ - 9 -- T- 

2r q 2r z„ 2r 2r 
Let u = 2r^lw^ to rewrite (V.23) more compactly as 


(V.24) 


u\"+{l + \-u)\'+ ^\ =0, (V.25) 

which has the form xy" + (/ + !- x)y' + py = 0 , of the differential equation for the 
associated Laguerre function [9], where I and p are real numbers. When p is a 

non-negative integer, i.e. /7 = 0,1,2,3,. Laguerre differential equation has as 

solutions the associated Laguerre polynomials L^(x). The first few polynomials 
are listed in Table 2. 



From equation (V.25), the Laguerre polynomials specify the constant ^ as 


2 


= P 


SO that 

w 


(V.26) 
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We can now use this result in the equation (V.20) in order to find the last 
unknown function P{t) 


0 

/ +1 // w ° 2 p ° i I i It 2i p 

- I P = - P = -1-1- 

Z,q w w" Z„q Z„q 


(V.27) 


We rewrite the last equation, using q(T) from (V.22) and = z„ + t^Iz„ , so that 


P - i Zg{\ + / + 2p )— - - - + 


2 2 2 2 ■ 

Z„ + T Z„ + T 


(V.28) 


Finally, integration from 0 to r gives 


( \ (+ 

P{t) = i (1 + / + 2p)arctarH^—j + ——-—j 


(V.29) 


At this point, we have found all the unknown functions P(t) , qir) and power 
series and we can put them into the trial solution (V.7). After rearranging 
terms, we can write down the solution to the parabolic wave equation of the (p, 1) 
mode. This solution is usually referred as the (p, /) Laguerre-Gaussian mode 
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V^rY 


V ^ y 




(V.30) 


with phase 


(r, i9,r) = -(1 + / + 2p) arctan 
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Z„ +T 
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(V.31) 


and the Laguerre polynomials listed in Table 2. 
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C. LAGUERRE-GAUSSIAN MODE ANALYSIS 

1. Intensity Plots 

Solving the parabolic wave equation in cylindrical coordinates, we find a 
family of modes, in equation (V.30), that are expressed in terms of the product of 
a Gaussian envelope, a Laguerre polynomial and a phase term. Each mode is of 
the {p, 1) order, where p is the radial mode number and / is the angular mode 
number. For this reason, the Laguerre-Gaussian modes are usually denoted as 
LG'p modes. 

Figure 15 shows the intensity profiles of various LG^ modes at r = 0, 

created from simulations using C code. In these plots there is no phase 
information displayed, since we are plotting the , x a* ^ quantity. Examining 

these plots, we are able to evaluate several features of the Laguerre-Gaussian 
beams. 

We see that the lowest order mode LG® is just the fundamental Gaussian 
mode, as expected. The other modes, with p = Q and />0, have a simple ring 
profile, since they have one central node. As the angular mode number I 
increases, the central, near-zero intensity region becomes larger, which is also 
true for modes with p>0. All modes with p>0 have a profile in the form of 
concentric rings. The radial mode number p indicates the number of nodes and 
therefore the number of the rings. Thus, the total number of dark nodes is p at 
r 0, in addition to the central node. These p + l nodes come from the zeros of 
the Laguerre polynomials. We also see that the width of the rings gets narrower 
for beams with higher I and same p. Finally, as the radial mode number p 
increases, the effective width of the beam gets bigger. 
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Figure 15. Laguerre-Gaussian modes LG^ ; Intensity plots (a^ , xa* ,) in the 
transverse plane x,y, at r = 0, where red is most intense, while 
blue is least intense. 


In order to see more clearly how an LG|, beam varies along the 

transverse plane, we plot in Figure 16 for two sample modes, the radial intensity 
distribution of the Laguerre-Gaussian beam across a diametric plane. 
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(a) (b) 

Figure 16. Radial intensity distribution for (a)LG 2 and (b)LG 3 beams. 


2. Plots of Real Part 

Plotting only the intensity of a Laguerre-Gaussian mode, we neglect 
phase. These effects turn out to play an important role in the evolution of a 
propagating beam. In Figure 17, we plot the same Laguerre-Gaussian modes as 
in Figure 15, but this time we include the phase term of equation (V.30) by 
plotting the real part squared. We could as well plot the real part, but in this case 
we would have to plot negative values of the optical field amplitude a. We 
immediately note that in all modes, except for the fundamental LG °, radial nodes 
appear. The number of nodes is proportional to the angular mode number /, 
since when r = 0, the phase (V.31) is just (r,i9,r) = /i9. 
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Figure 17. Laguerre-Gaussian modes LG[ ; plots of the real part squared, at 
r = 0. 
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D. PROPAGATION OF LAGUERRE-GAUSSIAN MODES 

We have seen previously in Figure 12 that a Gaussian beam, which is an 
LGI beam, retains its shape as it propagates. We find, by running similar 

simulations, that this is true for all higher-order modes, provided the laser beam 
consists of only one pure mode. 

1. Intensity Plots 

Figure 18 shows the results of computer simulations for the evolution of 
several Laguerre-Gaussian modes as they propagate from r = 0 to t = \. We 
show the intensity plots of the transverse cross section at the beginning and at 
the end of the propagation, and the cross section of the beam from r = 0 to r = 1. 
It is clear that as the beam evolves, its shape is unchanged. The central node is 
preserved, and even though the beam diverges, the multiple radial nodes persist. 
Nevertheless, we can see clearly that the peak intensities of the rings (antinodes) 
decrease as a function of the propagation distance, showing that the beam 
spreads reducing intensity. 
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Figure 18. Intensity plots of Laguerre-Gaussian modes as they propagate. 
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2 . 


Plots of Real Part 


In order to complete our study of the evolution of the Laguerre-Gaussian 
beams, we must examine the phase structure of the beam as it propagates. 
Simulations have shown that we have to take into account the well-known fact 
that light carries angular momentum [10]. For example, in the case of a circularly 
polarized beam, light caries an angular momentum that relates to the spin of 
individual photons. Furthermore, a light beam can also have orbital angular 
momentum, a feature that can be deduced classically from Maxwell’s equations. 
It has been shown that for a Laguerre-Gaussian beam, the orbital angular 
momentum is well defined and has a value of Ih per photon, where / is the 
angular mode number and h is the Planck’s constant [11]. 

The orbital angular momentum arises from the azimuthal dependence, 

particularly from the term of the solution to the wave equation (V.30). We 
will analyze this solution a bit more by examining a sample of Laguerre-Gaussian 
beams. We will plot the real part of the solution (V.30), so as to include the phase 
and see how it changes along the propagation distance. In the following figures 
we show the surface plots of a few Laguerre-Gaussian modes at the end of the 
propagation distance (where r=l). In this way, we can see how the phases 
evolve, compared to their initial shapes at r = 0, presented in Figure 17. For 
each mode we plot the real part and the real part squared, so that in the latter we 
avoid plotting negative values, as we did in Figure 17. 

Starting with the LGl mode, we put r = l in the solution (V.30) and we 
use Table 2 to find the associated Laguerre polynomial with p = 0. Then, the real 
part of the solution can be written explicitly as 
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Figure 19 shows two surface plots at time r = 1, for a LGl laser mode. In Figure 
19(a) we plot the real part, as described by equation (V.32), while in Figure 19(b) 
we plot the real part squared. We notice that a spiral pattern appears in the mode 
structure, due to the fact that the phase rotates as a function of r and the angle 
i9. We have to stress that the light does not follow a spiral path; it is just the 
phase that changes in a way that it describes such a spiral path. This path, in 
three dimensions is somewhat helical [10]. 



(a) (b) 

Figure 19. Plots of (a) real part and (b) real part squared of a LGl beam at the 
end of the propagation distance (r = 1), with = 1 and = 0.4. 


Other examples of the evolution of Laguerre-Gaussian beams, including 
phase change, are shown in Figure 20 and Figure 21, where we present, as in 
Figure 19, the surface plots of a LGf and a LG^ laser beam respectively. 



(a) (b) 

Figure 20. Plots of (a) real part and (b) real part squared of a LGl beam at the 
end of the propagation distance (r = 1), with a„ = 1 and z„ = 0.4. 
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The solution to the parabolic wave equation for the LG\ mode can be 
written in more detail as 



r (11 

— - -6arctan — +.9 ,(V.33) 

4 +1 K^o) \ 


giving the plots of Figure 21. We observe again the spiral patterns due to the 
phase evolution. We can see more clearly how varies with r in Figure 22, 

where we show the radial distribution of along various angles. 



(a) 


(b) 


Figure 21. Plots of (a) real part and (b) real part squared of a LG\ beam at the 


end of the propagation distance (r = 1), with = 1 and z„ = 0.4. 


Analyzing equation (V.33), we recognize the term e as the radial 
damping term. The behavior of the oscillating cosine term is easier to understand 
when we plot the set of points in r and & where it is maximum (Figure 23). We 
see that due to this term, the spiral pattern appears. These two terms, along with 
the term that comes from the Laguerre polynomial, give the beam profile of 
Figure 21(a). 


56 













(a) 


r 


(b) 



r 


(C) 

Figure 22. Radial distribution of a LG^ beam along (a) 3 = 0, (b) 3 = nll and 
(c) 3 = njA. 



Figure 23. Set of points in r and 3 where the cosine term of equation (V.33) 
is maximum. 
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VI. HIGH-ORDER MODES IN FREE ELECTRON LASER 

SIMULATIONS 


In previous chapters, we have studied the fundamental Gaussian optical 
mode as well as the higher-order modes in both rectangular and cylindrical 
coordinates and their propagation through a medium without any sources. The 
primary tool for our research was the computer code, which enabled us to solve 
the parabolic wave equation numerically. To analyze the operating conditions of 
an FEL inside the undulator, where the optical beam coincides with the electron 
beam, is more complex. There we must solve both the wave equation (111.52) and 
the pendulum equation (111.22) at the same time. At NPS, the FEL research group 
has developed computer codes that solve these two equations simultaneously, 
and they can be used to simulate various working parameters for a free electron 
laser system. These are powerful tools to gain deeper understanding of FEL 
physics and determine the importance of every operational parameter. Perhaps 
more importantly, it helps us to optimize FEL performance. 


A. ELECTRON BEAM TILT 

One must first understand the FEL under ideal operating conditions. But in 
order to cover more realistic situations, we must study the system performance 
under non-ideal conditions. Vibrations of a mobile system for instance, might 
cause the electron beam to enter the undulator off-axis or at an angle. We use 
the three-dimensional simulations (in x, y and r) developed in the FEL 
research group to investigate how sensitive the FEL is to off-axis tilting and 
shifting. In the simulations, we use dimensionless transverse coordinates 

X = XyJk/2L , y = yyJk/2L and the normalized time r = cr/L. Three-dimensional 

simulations were used to analyze the tolerance of the Jefferson Lab (Jlab) FEL to 
a tilted electron beam. The simulations discussed below, are divided into two 
parts. 
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1 . 


Gain 


In the first part of the Jlab simulations, we explore the gain response to the 
misaligned electron beam. In other words, we focus our attention on region 1 of 
Figure 24, which is a graph of a typical power build up of an FEL. In this region, 
we have just turned on the FEL so it has not yet reached steady-state. The gain 
is proportional to the slope of the power, in this case. 



Figure 24. Simulation output, showing the power P versus number of passes 
n in an FEL. In the first region, power is still building up while in 
region 2, FEL has reached its steady-state. 

We use the basic input parameters shown in Table 3, where we simulate 
weak fields (a^ = 71) and small number of passes n so as to be sure that we 
explore the power evolution in region 1. Also, the electron beam focus position is 
Tp = 0.5 , SO the electron beam is focused at the center of the undulator. We thus 

run our simulations by varying the tilt angle about the center of the undulator. 

We also normalize the tilt angle of the electron beam to the undulator 
length and the optical wavelength 


3 


yo 


3 

yo 


-- for the Jlab parameters. 

(426mrad) 


(VI.1) 
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Rayleigh Length 


0.42 

Initial Optical Field Amplitude 

a. 

0.001 

Current Density j 


20.0 

Number of Undulator Periods 

N 

30 

Number of Passes n 


30 

Electron Beam Focus Position 


0.5 

Optical Waist Position 


0.5 

Cavity Quality Factor 


17.0 


Table 3. The basic dimensionless parameters used in the Jefferson Lab FEL 
simulations for the Gain. 

We vary the normalized tilt angle from 0 to 4. To run the simulations, we 
enter the parameters of Table 3 and a range for the initial phase velocity . The 
optimum value of v; is determined in this way and we then look at the output of 
the run with this specific phase velocity v^. In Figure 25, we show a sample 
output for = 2.0, where the optimum initial phase velocity was found to be 
= 5.5. 

The dimensionless input parameters are shown in the upper-right box of 
Figure 25. On the upper-left graph, we plot the evolution of a slice of the optical 
beam at the outcoupling mirror, over the number of passes n , along with the 
color scale for the intensity plots of the optical field amplitude. Next to it, there is 
a 3-D representation of the beam shape after the last pass, again at the 
outcoupling mirror. Below the 3-D graph, we plot the cross-section of the optical 

beam |a(0,y, r)| (yellow) and the electron beam (red) inside the undulator, at the 
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last pass. Further down, there is a graph of the evolution of the phase velocity 
distribution over the number of passes /(i/,n), and next to it there is a phase- 
space plot of the electrons at the end of the undulator. Below the parameter 
table, there are two plots which show the development of the optical power P(n) 
and the gain G(n) . In the next row, there are graphs at both the outcoupling and 
the left mirror after the last pass, showing the optical beam transverse profile 
(yellow) compared to the theoretical Gaussian (purple), and the optical phase. 
Also shown in the same graphs is the intensity of the light (using the color scale 
at the top) and the extent of the mirrors (black). At the bottom, there is a list of 
the calculated results and in the lower right we plot the coefficients of the modes 
|c(m,n)| that comprise the beam at the outcoupling mirror after the last pass. The 

modal analysis is performed in terms of the Hermite-Gaussian modes. In the 
simulation results of Figure 25 we see that the output beam comprises basically 
from the fundamental mode (m = n = 0) but simulation predicts that there will be 
few higher-order modes in the beam structure as well. This is shown with the 
plotted squares (lighter blue color). Using this tool, we can look at the mode 
composition of any arbitrary wavefront that our simulations predict as the laser 
output. 

Figure 26 is a summary of the results showing gain G versus electron 
beam tilt angle We note that the FEL can operate pretty well with an 
imperfect electron beam injection into the undulator, since the gain falls to half of 
its peak value at dimensionless tilt angle of about « 3.5 , corresponding to an 
actual tilt of about 1.5 mrad. 
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nx=200, nt=50, np=30000, T^=8, seed=7 


Thu Feb 15 11:10:07 2007 


FEL imVEFRONT EVOLUTION 
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Figure 25. Output of a 3-D Simulation for Jefferson Lab FEL. 
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Electron Beam Tilt 


Figure 26. Summary of the simulation results for the gain versus increasing 
electron beam tilt angle 3^^. 


2. Extraction 

Having explored how the FEL responds to a tilted electron beam during 
the early stages of its operation, we can now explore the effects during its steady 
state phase, region 2 of the graph in Figure 24. So we choose strong initial 

optical fields and large number of passes n , as shown in the Table 4. 

As a measure of the performance of the FEL in this case, we estimate the 
extraction rj, defined as the fractional energy transferred from the electron beam 
to the optical field during one pass through the undulator 

extracted optical power 

. -■ (VI.2) 

initial electron beam power 
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Rayleigh Length 

0.42 

Initial Optical Field Amplitude 

10 

Current Density j 

20.0 

Number of Undulator Periods N 

30 

Number of Passes n 

500 

Electron Beam Focus Position Zp 

0.5 

Optical Waist Position z^ 

0.5 

Cavity Quality Factor 

17.0 


Table 4. The basic dimensionless parameters used in the Jefferson Lab FEL 
simulations for the extraction. 

Running the simulations as before, we first determine the value of the 
phase velocity v; for peak extraction rj, and then we look at the output of this 
particular run. One sample output for .9^^ = 3.0 is shown in Figure 27, where we 
can verify from the power evolution graph P(n), that the FEL is in its steady-state 
operation region. 
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Figure 27. Output of a 3-D Simulation for Jefferson Lab FEL, exploring 
extraction. 
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Summarizing the results for many tilt angles we note with some 
surprise that the FEL exhibits good tolerance; the extraction actually increases 
with increasing angles up to a point, as shown in Figure 28. This may be due 

to “overbunching” of the electrons when there is no tilt, reducing the extraction in 
that case. The simulation predicts extraction greater than 1% for normalized 
beam tilt corresponding to an actual tilt about l.Smrad. 



Electron Beam Tilt 


Figure 28. Summary of the simulation results for the extraction ?] versus 
increasing electron beam tilt angle . 
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VII. CONCLUSION 


In this thesis, after reviewing the basic FEL and optical theory, we found 
solutions to the parabolic wave equation that suggest the existence of higher- 
order optical modes, working with either rectangular or cylindrical coordinates. 
Motivated by the symmetry of the FEL system, we focused our attention on the 
solutions derived using cylindrical coordinates, which employed Laguerre 
polynomials. Analyzing the Laguerre-Gaussian modes, we found interesting 
results for the beam shape and for the beam evolution as well. The modal 
analysis showed us that we have to consider the physical attributes of light, in 
order to explain the unexpected behavior of the laser beam after its propagation. 
Particularly, we noticed spiral patterns in the beam, attributed to the orbital 
angular momentum of light. Finally, we demonstrated the appearance of higher- 
order modes in an FEL output, while investigating the tolerance of the Jefferson 
Lab FEL to electron beam tilt. 

Currently, the NPS FEL research group is analyzing the output wavefronts 
of computer simulations in terms of Hermite-Gaussian modes. We suggest that a 
future thesis work would develop a function that analyzes an arbitrary wavefront 
from an FEL simulation in terms of the Laguerre-Gaussian modes. The use of 
this tool could make the optical beam analysis of any simulated FEL output more 
complete, helping us understand better the effects of the investigating operating 
parameters. 
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